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Shocked Stellar Wind or Collimated Outflows? * 
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Abstract. The Eskimo Nebula (NGC 2392) is a double-shell planetary nebula (PN) known for the exceptionally large expan- 
sion velocity of its inner shell, ~90 km s -1 , and the existence of a fast bipolar outflow with a line-of-sight expansion velocity 
approaching 200 km s . We have obtained XMM-Newton observations of the Eskimo and detected diffuse X-ray emission 
within its inner shell. The X-ray spectra suggest thin plasma emission with a temperature of ~ 2 x 10 6 K and an X-ray lumi- 
nosity of Lx = (2.6 ± 1.0) x 10 3i 1 150 pc) 2 ergs s" 1 , where d is the distance in parsecs. The diffuse X-ray emission shows 
noticeably different spatial distributions between the 0.2-0.65 keV and 0.65-2.0 keV bands. High-resolution X-ray images of 
the Eskimo are needed to determine whether its diffuse X-ray emission originates from shocked fast wind or bipolar outflows. 
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1. Introduction 

Planetary nebulae (PNe) consist of the stellar material ejected 
by low- and intermediate-mass stars (1-8 M ) at the end of the 
asymptotic giant branch (AGB) phase. As such a star evolves 
off the AGB, the copious mass-loss strips off the stellar enve- 
lope and exposes the hot stellar core. A PN emerges when the 
stellar UV radiation ionizes the ejected stellar material, causing 
it to emit in the optical. 

PNe are expected to be diffuse X-ray sources. The central 
stars of PNe possess fast st ellar winds with termina l veloc - 
ities of 1000-4000 km s 1 dCerruti-Sola & Perinottol Il985t> . 
while fast collimated outflows with velocities up to 1000 
km s _1 are also ob served in PNe and proto-PNe, e.g., MvCn 18 
and Hen 3-1475 (iBobrowskv et all M99.4 iRieraetall 1 199.4 
lO' Connor et all l2000t) . The interactions of the fast stellar 
wind and/or collimated outflows with nebular material produce 
shocked gas that is hot enough to emit in X-rays. 

In the interacting-ste l lar-win ds model of PNe 
jKwok. Purton. & F itzgerald. Il978h . the fast wind ema- 
nating from the central star sweeps up the slow AGB wind to 
form a sharp nebular shell. The interior structure of a PN would 
be similar to that of a wind-blown bubble (e.g.. IWeaver et all 



I19771 The central cavity of a PN is expected to be filled 
with shocked fast wind at 10 7 — 10 8 K, but this hot gas is too 
tenuous to produce appreciable X-ray emission. Dynamic or 
evaporative mixing of cool nebular material into the hot gas 
at their interface produces optimal conditions for soft X-ray 
emission, which will show a limb-brightened morphology 
withi n the nebular shell, as observed in NGC 6543 dChu et all 
l200ll) . 

Fast collimated outflows or jets may also produce hot X- 
ray-emitting gas. When outflows with velocities >300 km s 
initially impinge on the AGB wind, bow sho cks and X- 
rav e mission can be produced (e.g., Hen 3-1475. ISahai et all 
120031) . The prolonged action of collimated outflows may bore 
through the AGB wind and form extended cavities filled b 
shock ed hot gas that emits X-rays (e.g., Mz 3. iKastner et 
2003). 
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ESA Member States and NASA. 



Diffuse X-ray emission fro m hot gas in PNe wa s hinte d 
by ROSAT observations llGuerrero. Chu. & Gruendl l2000h . 
but was unambiguously resolved only by Chandra and 
XMM-Newton observations. Besides the aforementioned 
Hen 3-1475, Mz3, and NGC 6543, diffuse X-ray emission 
has been reported in only 4 other PNe: A 30, BD+30°3639 
NGC 7009, and NGC 702 7 (IChn Chang. & Con wavi " 
iGuerrero. Gruendl. & Chul 120021 IKastner et all 



1997 



2000; 



Kastner. Vrtilek. & Soker, 2001). It is imperative to detect dif- 
fuse X-ray emission from a large number of PNe to investigate 
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Fig. 1. XMM-Newton and HST Ha images of the Eskimo Nebula. Panel (a) displays the XMM-Newton EPIC raw image in the 
0.2-2.0 keV band with pixel size 1'.'5. Panel (b) shows the HST Ha image overplotted by the 25%, 50%, 75%, and 95% X-ray 
contours extracted from the adaptively smoothed 0.2-2.0 keV band image described in the text. Contours at fainter levels are not 
plotted, as they show mainly the extended wings of XMM-Newton's PSF and the smear caused by the adaptive smoothing. Panels 
(c) and (d) display XMM-Newton EPIC raw images in the 0.2-0.65 and 0.65-2.0 keV energy bands, respectively. The pixel size 
in these images is 1" 5. The arrows in panel ( d) mark th e location of the fast collimated outflow as derived from high-dispersion 
echelle spectroscopy (Gieseking, Becker, & Solf, 1985, Maestro et al., in preparation). 



whether shocked fast stellar wind, collimated outflows or both 
are responsible for the X-ray emission from PNe. 

The Eskimo Nebula (NGC2392) offers an excellent op- 
portunity to assess the effects of collimated outflows in a PN. 
The Eskimo is a double-shell PN: its bright elliptical inner 
shell has two blister-like protrusions on the south rim, and its 
round outer shell contains a set of complex low-ionization fea- 
tures JO'Dell etall l2002h . In addition, the Eskimo has a fast 
collimated outflow, detected in high-dispersion spectra of op- 
tical nebular lines; the outflow, with line-of-sight expansion 



velocities approaching 200 
ing with nebular material 


km s \ is dynamically interact- 


Gieseking. Becker. & Solf, 


1985; 


lO'Dell & Balllll985HO'Dell. Weiner. & Chulll990. Maestro et 



al., in preparation). We have obtained XMM-Newton observa- 
tions of the Eskimo and here we report the detection of diffuse 
X-ray emission from this nebula. 



2. Observations 

The Eskimo Nebula was observed with the XMM-Newton 
Observatory in Revolution 790 on 2004 April 2 using the 
EPIC/MOS1, EPIC/MOS2, and EPIC/pn CCD cameras. The 
two EPIC/MOS cameras were operated in the Full-Frame 
Mode for a total exposure time of 17.5 ks, while the EPIC/pn 
camera was operated in the Extended Full-Frame Mode for 
a total exposure time of 12.5 ks. For all observations, the 
Medium filter was used. The XMM-Newton pipeline prod- 
ucts were processed using the XMM-Newton Science Analysis 
Software (SAS version 6.0.0) and the calibration files from the 
Calibration Access Layer available on 2004 June 10. 

The event files were screened to eliminate events due to 
charged particles or associated with periods of high back- 
ground. For the EPIC/MOS observations, only events with 
CCD patterns 0-12 (similar to ASCA grades 0—4) were se- 
lected; for the EPIC/pn observation, only events with CCD 
pattern (single pixel events) were selected. Time intervals 
with high background (i.e. count rates > 0.3 cnts s for 



the EPIC/MOS or > 1.4 cnts s" 1 for the EPIC/pn in the 
background-dominated 10-12 keV energy range) were dis- 
carded. The resulting exposure times are 17.3 ks, 17.1 ks, and 
1 1.6 ks for the EPIC/MOS 1, EPIC/MOS2, and EPIC/pn obser- 
vations, respectively. 

3. Results 

The XMM-Newton EPIC/MOS 1, EPIC/MOS2, and EPIC/pn 
observations of the Eskimo detect within the nebula a total of 
180 + 15, 175 + 15, and 620 + 30 cnts, respectively. In order 
to construct an X-ray image of the highest signal-to-noise ra- 
tio, we merged together the event files of the three EPIC ob- 
servations and extracted raw EPIC images in the 0.2-2.0, 0.2- 
0.65, and 0.65-2.0 keV bands with a pixel size of 1'.'5 (Fig. la, 
lc and Id). The raw EPIC image in the 0.2-2.0 keV band is 
then adaptively smoothed using Gaussian profiles with FWHM 
ranging from 4" to 6". The contour map of this smoothed im- 
age is overplotted on a Hubble Space Telescope (HST) Wide 
Field Planetary Camera 2 (WFPC2) Ha image (Fig. lb). The 
alignment of X-ray and optical images is fine-tuned using 
HD 59087, a star ~100" north of the Eskimo, which is detected 
both in the XMM-Newton EPIC and HST WFPC2 images. 

The X-ray emission from the Eskimo is clearly extended. 
Its distribution is elongated along PA~25° (Fig. la), and the 
25% contour of the smoothed EPIC image follows closely the 
outline of the the inner shell in the HST Ha image (Fig. lb). 
The image in the 0.2-0.65 keV band shows a similar spatial 
distribution, but in the 0.65-2.0 keV band the emission is elon- 
gated along PA~70°, i.e., roughly aligned with the fast bipolar 
outflow detected at PAs of 50°-80° and 230°-260 °, as illus- 
trated in Fig. Id JGieseking. Becker. & Sold fl985l Maestro et 
al., in preparation). These comparisons suggest that the diffuse 
X-ray emission from the Eskimo is mostly confined within its 
inner shell, but some of the harder X-ray emission in the 0.65- 
2.0 keV band may be produced by the interaction of the fast 
bipolar outflow with nebular material. 
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Fig. 2. (top) XMM-Newton EPIC/pn background-subtracted 
spectrum of the Eskimo overplotted with the best-fit MEKAL 
model. For plotting purposes, both the spectrum and model are 
binned in 30 eV wide energy bins. ( inset) x 2 grid plot as a func- 
tion of kT and A^h of the spectral fit to the EPIC/pn spectrum. 
The contours represent 68%, 90%, and 99% confidence levels. 
(bottom) Relative residuals of the best-fit model to the EPIC/pn 
spectrum. 



To analyze the properties of the diffuse X-ray emission 
from the Eskimo, we extracted a spectrum from the EPIC/pn 
event file using a circular source aperture of 20" radius, large 
enough to include all X-ray emission from the nebula. The 
background level was assessed from an annular region exte- 
rior to the source aperture and 9 times larger in area. The 
background-subtracted EPIC/pn spectrum of the Eskimo is 
shown in Fig. 2. This spectrum shows a broad peak between 
0.4 and 0.6 keV, a plateau below 0.4 keV, and a sharp decline 
above 0.6 keV. Low-intensity emission peaks are discernible at 
~0.9 keV, ~1.1 keV, and -1.8 keV. The overall spectral shape 
indicates thermal plasma emission, and the peaks correspond 
to the He-like triplets of N vi at -0.43 keV, O vn at -0.57 keV, 
Ne ix at -0.92 keV, Na x at - 1 . 1 3 keV, and Si xm at - 1 .86 keV. 
We have also extracted EPIC/MOS spectra; they show similar 
spectral features at lower signal-to-noise ratios, so we will con- 
centrate only on the EPIC/pn spectrum in the analysis below. 

The observed spectrum can be modeled to determine the 
physical conditions of the X-ray-emitting gas and the amount 
of foreground absorption. We have adopted the MEKAL 
optically-thin plasma emission model jjKaastra&MeweLll993t 
iLiedhal. O sterheld. & Golds tein! 1 19951) and nebular abun- 
dances, although we note that th ese may diffe r from those 
of the X-ray-emitting gas (e.g., iManess et all l2003h . The 
abundances of He, C, N, O, Ne, Ar, and S relative to 
hydrogen by number, as determined from optical and UV 
spectra of the Eskimo, are 0.076, 2.2 x 10~ 4 , 1.1 x 10~ 4 , 
2.8 x KT 4 , 6 4 x 1(T 5 , 1.4 x 10~ 6 , and 4.3 x 10^ , respec- 
tively llBarkerl Il99lt iHenrv. Kwitter & BatesL l2000l) . For el- 
ements that do not have available nebular abundances, we 



adopt solar abundances JGrevesse & Sauval Il998h . We have 
also assumed solar abundances for the foreground interstel- 
lar a bsorptio n, and adopted absorption cross-sections from 
iBalucinska-Chu rch & M cCammonl dl992h . 

The spectral fit is carried out by folding the absorbed 
MEKAL model spectrum through the EPIC/pn response ma- 
trix, and comparing the modeled spectrum to the observed 
EPIC/pn spectra in the 0.2-2.5 keV energy range using the 
X 2 statistics. The best-fit model with nebular abundances has 
a large reduced x 2 , -2.5, with excessive positive residuals at 
-0.4 keV and -0.9 keV and negative residuals at -0.6 keV. 
These energies are coincident with the N vi, Ne ix, and O vn 
lines, respectively; thus the residuals suggest that the N/O and 
Ne/O abundance ratios are higher than those of the adopted 
nebular values. We have therefore made spectral fits that al- 
lowed the abundances of N and Ne to vary, and subsequently 
the reduced^ 2 is improved to -1.5. 

The best-fit model, overplotted on the EPIC/pn spectrum 
in Fig. 2, has a plasma temperature of T = (2.0^ 3) x 10 6 
K (or kT = 0.175^q[^ keV), an absorption column density 
Nh = (8^i 2 ) x 10 20 crrT 2 , and a volume emission measure of 
(1.2+0.3)xl0 5 V/1150pc) 2 cm" 3 , where d i s the distance in 
parse cs and has been reported to be 1 150 pc dPerrvman et all 
1997). The quality of the spectral fits is illustrated by the plot 
of the reduced x 1 of the fits as a function of kT and Nh shown 
in the inset of Fig. 2. A dopting a typical gas-to-dust ratio 
(Bohlin. Sava ge. & Drake! Il978l) . the best-fit absorption col- 
umn density corresponds to a logarithmic extinction at the H/3 
line of ciy3 = 0.2^-j, in agreement wi th that determined from 
the B aimer decrement (Barker, 1991). The N and Ne abun- 
dances are enhanced with respect to the nebular abundances by 
3.5+1.2 and 3.0+1.0, respectively, thus resulting in N/O-1.4 
and Ne/O~0.7. The residuals of the best-fit model may also 
indicate enhanced Na and Si abundances, but the number of 
counts is too small to warrant spectral fits allowing these abun- 
dances to vary. 

The observed (absorbed) X-ray flux in the 0.2-2.5 keV en- 
ergy range is (6 ± 1) x 10~ 14 ergs crrr 2 s _1 , and the intrinsic 
(unabsorbed) X-ray flux is (1.6 + 0.6) x 10~ 13 ergs crrr 2 s -1 . 
The X-ray luminosity in the same energy range is Lx = (2.6 + 
1.0) x 10 31 (<i/1150pc) 2 ergs s" 1 . 

4. Discussion 

The sharp rim of the Eskimo's inner shell suggests compression 
by supersonic shocks. Indeed, the inner shell of the Eskimo, ex- 
panding at -90 km s _1 into the much slower ou ter shell, has an 
expansion velocity notably high among PNe JO'Dell & BallL 
1985). To assess whether the hot gas in the central cavity is 
responsible for driving such a fast expansion, we derive the 
thermal pressure of the X-ray-emitting gas and compare it to 
those of the surrounding nebular shells. Assuming a prolate 
ellipsoidal central cavity, the volume occupied by the X-ray- 
emitting gas is - 1.1 x 10 51 (<i/l 150 pc) 3 (e/0.5) cm 3 , where e 
is the filling factor and may be -0.5. From this volume and the 
aforementioned volume emission measure, we derive an rms 
electron density of « e -35 (d/l 150pcT I/2 (e/0.5)" I/2 cm" 3 . 
This rms electron density and the plasma temperature of 2.0 x 



4 



M. A. Guerrero et al.: Hot Gas in the Eskimo Nebula 



10 6 K imply that the thermal pressure of the X-ray-emitting 
gas, P th ~ 1.9 x n e kT, is ~ 2 x 10~ 8 (d/l 150pc)- 1/2 (e/0.5r 1/2 
dynes cm 2 . The inner and outer shells of the Eskimo have den- 
sities of 2500 and 900 cm " 3 , respectively, and a temperature 
of 14,500 K dBarkerlll99lT) . The thermal pressure of the inner 
and outer shells are ~lxl0~ 8 and ~3xl0~ 9 dynes cm -2 , respec- 
tively. The X-ray-emitting gas has slightly higher thermal pres- 
sure than the inner shell, and much higher than the outer shell, 
thus the thermal pressure of the hot gas in the Eskimo's central 
cavity drives the expansion of its inner shell into the outer shell. 

The temperature of the X-ray-emitting gas in the Eskimo, 
2.0xl0 6 K, is similar to that found in other elliptical PNe; 
however, the terminal vel ocity of the Eskimo's fast w ind 
is low , only ~400 km s" 1 JPauldrach. Hoffmann. & Mendezl 
120031) . For an adiabatic shock, the expected temperature of the 
shocked wind would be ~1.9xl 6 K for such a velocity, barely 
reaching the temperature indicated by the X-ray spectrum, in 
sharp contrast to other PNe whose shocked stellar win ds are 
10-100 times hotter than th e hot gas detected (e.g.. lChu et all 
l200ll iKastner et all l2000l) . If the hot gas in the Eskimo in- 
deed originates from the shocked fast stellar wind, the low ex- 
pected post-shock temperature implies that no significant mix- 
ing of nebular material has taken place. The mass of the hot 
gas, ~3xl0~ 5 (d/l 150pc) 5/2 (e/0.5) 1/2 M , can be supplied by 
the fast stellar wind in ~180 yr at a constant mass loss rate of 
1.8X10" 8 M yr 1 dPauldrach. Hoffmann. & Mendezl 12003b . 

If the X-ray emission from the Eskimo does not originate 
from a shocked fast stellar wind, then there are two other pos- 
sibilities. First, the X-ray emission can be partially attributed 
to the central star as in the case of NGC 6543 o r to a late-type 
binary companion as in the case of NGC 7293 JGuerrero et all 
1200 lb . The angular resolution of XMM-Newton, however, is in- 
sufficient to resolve a point source from the diffuse emission 
from the Eskimo. Second, the diffuse X-ray emission can be 
produced by the dynamical interaction of the fast bipolar out- 
flow with the inner shell. This collision will produce shock- 
excited gas along the direction of the outflow which is de- 
tected as the harder X-ray emission in the 0.65-2.0 keV band, 
and al so as shock-excited r Fe nl 1 .26 and 1 .64 pm line emis- 
sions llHollenbach & McKeei Il989t iHora. Latter. & Deutschl 
119991) . Similar situation is observed in BD+30°3639, whose 
diffuse X-ray emission shows an asymmetric spatial distribu- 
tion in hard energies an d reveals spec t ral evi dence of enhanced 
Ne/O abundance ratio ( Kastner et al. l | 2002l) . As suggested for 
BD+30°3639 JSoker & Kastner! l2003h . die hottest gas in the 
Eskimo may have been produced by a fast, ~500 km s _1 , colli- 
mated post-AGB wind that is no longer present. Alternatively, 
the observed 200 km s outflow may have a large inclination 
with respect to the line-of-sight, so that its real expansion ve- 
locity is large enough to power the hottest gas in the Eskimo. 

Our XMM-Newton observations of the Eskimo Nebula have 
detected diffuse X-ray emission from hot gas within its central 
cavity. The high thermal pressure of this hot gas is responsible 
for the high expansion velocity of its innermost shell. The ori- 
gin of this hot gas, however, is uncertain, and may consist of 
several components: shocked fast stellar wind, shocks associ- 
ated with the fast bipolar outflow, and emission from its cen- 
tral star or a binary companion. X-ray observations at the high- 



est spatial resolution afforded by Chandra are needed to deter- 
mine accurately the origin of X-ray emission from the Eskimo 
Nebula, a very interesting and complex PN. 
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